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Being able to reliably track perturbations across bounces and turnarounds in cyclic and bouncing 
cosmology lies at the heart of being able to compare the predictions of these models with the 
Cosmic Microwave Background observations. This has been a challenging task due to the unknown 
nature of the physics involved during the bounce as well as the technical challenge of matching 
perturbations precisely between the expansion and contraction phases. In this paper, we will present 
general techniques (analytical and numerical) that can be applied to understand the physics of 
the fluctuations, especially those with “long” wavelengths, and test its validity in some simple 
bouncing/cyclic toy models where the physics is well understood. We will then apply our techniques 
to more interesting cosmological models such as the bounce inflation and cyclic inflation. 


I. INTRODUCTION 

It has long been appreciated that Standard Big Bang 
(BB) cosmology cannot be the complete theory of the 
early universe since it is plagued by an initial singu¬ 
larity [l|. The singularity theorems have also been ex¬ 
tended to apply to scalar field-driven inflationary cos¬ 
mology [2|, Q. It is possible that the initial singularity is 
resolved by quantum gravity where the space-time con¬ 
tinuum description is replaced by a more fundamental de¬ 
scription. Another possibility, however, is that an “effec¬ 
tive” space-time description is possible even at Planckian 
energy densities and corrections to the gravitational La- 
grangian and/or the description of matter stress-energy 
tensor leads to a bouncing universe , for a recent dis¬ 
cussion see {?!■ Traditionally, bouncing cosmology has 
been a popular alternative to the inflationary paradigm 
as during the contraction phase it is relatively easy to 
arrange for the Hubble radius H~ 1 (t ) to decrease faster 
than the physical wavelength of fluctuations. Hence, it is 
conceivable that processes acting in the contracting phase 
can lead to a non-inflationary mechanism for generating 
superHubble fluctuations that can explain the origin of 
structure in the universe. Proposals which envisage such 
a mechanism include the Pre-Big-Bang [|J, Ekpyrotic Q, 
Matter bounce mm and Hagedorn bounce fl~2l - [~Ml ) . 

More recently, it has been recognized that bounce and 
inflation need not be mutually exclusive leading to var¬ 
ious Bounce-inflation scenarios jTHi - fTof which opens up 
new ways to construct phenomenologically viable non¬ 
singular and geodesically complete models. For instance, 
in [151] it was recognized that since a bounce necessarily 
implies a super-inflationary phase preceding ordinary in¬ 
flation, the spectrum of modes exiting during the bounce 
would have a slight blue tilt as opposed to the conven¬ 


1 Another nonsingular way of resolving the BB singularity is to 
have an emergent universe 11 type of scenario where in the 
infinite past our universe approaches a Minkowski space-time, 

see also IQ]. 


tional red tilt that is generated during the inflationary 
phase. This can potentially address the increase in power 
that has been observed in the low (Z < 40) multipoles in 
CMB|H[2l|. 

Another class of bounce inflation scenarios have been 
considered in jH . [22I ] where the goal was to find a new 
mechanism that can yield a near-scale-invariant primor¬ 
dial spectrum during the contraction phase in the “con¬ 
stant” mode (C c mode) involving string thermodynam¬ 
ics. To be consistent with the observed amplitude of fluc¬ 
tuations in CMB, typically the fluctuating modes in these 
scenarios must exit the Hubble radius at very high tem¬ 
peratures ~ 10 15 GeV. This also means that the physi¬ 
cal wavelengths of these fluctuations are very small and 
needs to be stretched by about 60 efolds to be cosrno- 
logically relevant today. It then becomes imperative to 
invoke an inflationary phase following the contraction 
and bounce phases to make the models viable. How¬ 
ever, an advantage of such a scenario is that the infla¬ 
tionary regime need not be near-exponential, any accel¬ 
erating phase, such as power-law inflation, would be suf¬ 
ficient. This may aid in building phenomenological mod¬ 
els from fundamental theories such as string theory. It 
is worth pointing out that in other bouncing scenarios, 
such as in ekpyrotic models, the mechanism for produc¬ 
ing scale-invariant perturbations proceeds via a differ¬ 
ent route: One typically produces a near scale-invariant 
power-spectrum in the “growing” mode (G mode) during 
contraction. During the exit the amplitudes are “small” 
and the wavelengths are “large”, but in the subsequent 
contraction phase the fluctuations are amplified to the 
observed ~ 10 _1 ° level. An additional mechanism is then 
introduced to transfer these fluctuations to the adiabatic 
constant mode. 

A key aspect to assess the viability of all these different 
models that emerges from the above discussion is about 
how the fluctuations in the “constant” and “growing” 
modes in the contraction phase is transferred to the con¬ 
stant mode (G e mode) and “decaying” mode (D mode) in 
the expanding phase. In particular, what matters is the 
spectrum in the constant mode of the expanding branch 
because that is what we see in CMB. For instance, for 
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the success of the hagedorn bounce model [HJ , we would 
want the near scale-invariant C c -mode spectrum to be 
transferred to the C e mode. On the other hand, it would 
be advantageous for the ekpyrotic or matter bounce mod¬ 
els if the G mode spectrum could be transferred to the 
C e mode. This paper is primarily devoted to developing 
general strategies and techniques, both analytical and 
numerical, to study how the fluctuations are transmitted 
across the bounces. 

Additionally we will also discuss evolution of perturba¬ 
tions in cyclic cosmological models which are extensions 
of the bouncing models that not only entail a bounce 
phase but also turnarounds (transitions from contrac¬ 
tions to expansions). In cyclic models such as the ekpy¬ 
rotic cosmology, it is not particularly important to track 
the perturbations across the turnarounds as the fluctua¬ 
tions during the turnaround are too small to be relevant 
for the next cycle, and “new” perturbations are generated 
at each cycle. In other models [23l - (28l | however, perturba¬ 
tions that we see in CMB today are generated over many 
many short cycles and therefore it is important to be 
able to track the fluctuations through the entire course 
of the cycle including the bounce and the turnaround. 
Therefore, in this paper our aim is to be able to track 
fluctuations across both turnarounds and bounces, and 
indeed over many cycles. We have succeeded in devel¬ 
oping robust analytical methods to track the two usual 
modes of superHubble fluctuations, and have tested our 
methods using numerical computations. In the process 
we also found certain general trends that subHubble fluc¬ 
tuations seem to follow, but the robustness of these latter 
trends will need to be investigated further in future. 

Before delving into the details of our framework, let 
us briefly discuss some of the challenges one faces in try¬ 
ing to understand the evolution of fluctuations in bounc¬ 
ing/cyclic cosmologies. While a vast body of literature 
already exists on the subject of perturbations around 
bouncing backgrounds, it is fair to say that we have rela¬ 
tively few general and conclusive results. The biggest 
hurdle to accomplishing this task in a comprehensive 
manner is simply the fact that one doesn’t yet have a 
theory of Quantum Gravity (QG), or even an effective de¬ 
scription of it that one can trust near the bounce which 
is expected to occur at the Planck scale 2 . Often the 
bounce mechanism is mediated by unphysical and unsta¬ 
ble ghost degrees of freedom, and one of the reasons dif¬ 
ferent groups have found different results when it comes 
to tracking the perturbations could be because the un¬ 
derlying physical models suffer from different inconsis¬ 
tencies/instabilities and therefore cannot completely be 
trusted. 

To by pass our ignorance of the precise bounce mecha¬ 
nism often, the evolution of the background and pertur¬ 
bations are modeled by a contracting GR phase that is 


2 For progress in this direction see [29(l37ll , for instance. 


matched to an expanding phase of Standard Big Bang 
cosmology through an instantaneous and often singular 
transition along a space-like hypersurface. For instance, 
one proposal has been to use the analog of the Israel 
matching conditions (matching conditions [38| which de¬ 
scribe the merger of two solutions of the Einstein equa¬ 
tions along a time-like hypersurface). These equations 
were discussed in [H, Ho[ • As pointed out by , dif¬ 

ferent matching prescriptions, however, can lead to dif¬ 
ferent results thereby questioning the reliability of these 
prescriptions. 

In this paper, our strategy will be to first focus on re¬ 
sults that can be applied to a large class of bouncing mod¬ 
els. We will see that in a separate universe approach |48j 
(see also [H, [50] it is possible to track two modes of su¬ 
perHubble fluctuations as long as we know the bouncing 
background, the precise knowledge of the new physics 
that makes a bounce possible will not be needed as long 
as some basic symmetries are respected. We will test 
our results by looking at toy bouncing and cyclic models 
that can be constructed within the General Relativistic 
framework, and therefore does not suffer from any ambi¬ 
guities coming from unknown new physics. We will then 
apply our framework to understand the evolution of su¬ 
perHubble perturbations in more realistic cosmological 
models, such as the emergent cyclic inflation 6j, [5ll456| 
and bounce inflation jl5Bl9| . In the process we will also 
be able to shed some light on how subHubble fluctua¬ 
tions evolve in cyclic cosmologies using numerical com¬ 
putations. 

The paper is organized as follows: In section El we 
will develop the analytical formalism to track super 
Hubble fluctuations exactly given the knowledge of the 
background evolution. In section IIII1 we will consider 
simple General Relativistic toy models which exhibits 
turnarounds and bounces and will obtain the evolution 
of perturbations numerically. We will then test the an¬ 
alytical results obtained in section El by comparing the 
numerical results in the long wavelength limit. Next, in 
section EY1 we will apply our formalism to study simple 
QG-mediated bouncing cosmologies and, in particular, 
demonstrate how one can obtain the transfer matrix re¬ 
lating the decaying and constant mode of the expanding 
branch to the constant and growing mode of the contract¬ 
ing branch. In section El we will use our results to inves¬ 
tigate the stability and viability of the emergent-cyclic- 
inflation models. Subsequently, in section rvTl we will look 
at the perturbations in bounce inflation models which 
will also show us ways to go beyond the formalism devel¬ 
oped as well as provide insights into the evolution of the 
subHubble fluctuations. We will conclude in section eh 
by summarizing our main results and pointing out future 
research directions. Our appendix section fYIIIl contains 
an illustrative example on how it is possible to derive 
the perturbation equations for superHubble fluctuations 
from the knowledge of the background evolution equa¬ 
tions. 
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II. ANALYTICALLY TRACKING 
SUPER-HUBBLE FLUCTUATIONS 


Our aim in this section is to provide a general algo¬ 
rithm to analytically track certain superHubble modes. 
The key physics ingredient which allows us to provide 
such an algorithm is the realization that for superHub¬ 
ble fluctuations one can ignore the spatial fluctuations in 
comparison with the time variation. In other words, in 
the evolution equations of the superHubble fluctuations 
one only need to keep track of the time derivatives, but 
this is precisely what is done to arrive at the cosmologi¬ 
cal background solutions. Thus, it should be possible to 
obtain the evolution equations of the superHubble fluc¬ 
tuations even if we only know the background evolution 
equation. In fact, if we already know the complete set 
of background solutions, then we can obtain the super¬ 
Hubble modes directly from them, without ever having to 
even know what the superHubble perturbation equations 
are! 


This is really a variant of the separate universe ap¬ 
proach. In the context of bouncing cosmology this ap¬ 
proach was considered in 57], but we provide a more 
general and comprehensive formalism with illustrative 
examples which also lets us verify the procedure with 
numerical precision. Our approach also provides several 
physical insights that could be helpful for model building- 
purposes. 


A. Super Hubble Modes from Background 
Evolution 


In the longitudinal gauge , in the absence of anisotropic 
stress, the perturbed metric is given by 


db(r) = a(r, {A,; = 0}), is now given by 

a(f) = a b { f)[l + ^2 Awhere 

, aoAi A . =0 

(IL3) 

Clearly, /)’s are nothing but the different superHubble 
modes with A^’s the corresponding amplitudes. However, 
the fluctuations have been obtained in the so called syn¬ 
chronous gauge 

ds 2 = a(f) 2 (1 + 2 ip) [— df 2 + dx 2 ] , (II.4) 

with ip = ^2 Kh , (II.5) 

i 

where for notational simplicity we have now dropped the 
subscript “b”. One needs to perform appropriate gauge 
transformations to obtain the more commonly discussed 
Newtonian potential, <h, in the longitudinal gauge (III. II) . 

To find the relation between the two gauges we perform 
the following transformations: 

ds 2 = a(r) 2 [(1 + 2<h)dr 2 — (1 — 2$)dx 2 ] 

« a(r) 2 (l—2<f>) [(1 + 4$)dr 2 - dx 2 ] = a 2 (l-24>) [df 2 -dx 2 ] 
where r and f are related as 

^ » 1 + 2$ (II.6) 

dr 

Although the above metric is now in the same form as 
(III. 41) , we see that there is an ambiguity in distinguishing 
the background scale factor from the perturbations. This 
has to be fixed by requiring form invariance of the scale 
factor. If £ relates the two time coordinates 


ds 2 = a(r) 2 [(1 + 2$)dr 2 - (1 - 2$)dx 2 ] , (II.l) 

where the relativistic potential $(x, r) is the field describ¬ 
ing the (scalar) fluctuations. The idea is to be able to find 
the superHubble fluctuations in $ from the knowledge of 
cosmological background solutions, 


f = r + £ => £' = 2$ (II.7) 

then 

a(r) = a(f) — £a'( f) = a(f)(l — £77) (II.8) 

This gives us the relation between $ and ip: 


ds 2 = a(T) 2 [—dr 2 + dx 2 ] , (II.2) 


ip = —$ - £,H (II.9) 


that satisfy the appropriate cosmological equations. 
Since the conformal times associated with the metrics 
(III. II) and (III. 31) need not be the same, we have intro¬ 
duced a “bar” to distinguish the two. Suppose now that 
the set of all cosmological solutions are known and are 
characterized by the parameter set {Aj}. These param¬ 
eters can be thought of as integration constants that 
arise when we solve the cosmological differential equa¬ 
tions. Let us further suppose that we are interested in 
finding perturbations around a specific background solu¬ 
tion chosen by convention to correspond to A,; = 0, V i. 
A small perturbation around the background metric, 


Thus, once we know ip from the cosmological solutions 
via (III.3|) and (III. 51) . we can obtain the gauge parameter 
from the equation, 

= (ii. io) 

derived using (III.91) and (III. 71) . Subsequently, we can 
compute the Bardeen potential $ from (III. 71) . In prac¬ 
tice, it may not be possible to know all the cosmological 
solutions, but rather a subset of them. This can then be 
used to obtain a subset of the superHubble modes, as we 
will explicitly illustrate in the next subsection. 
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B. Super-hubble modes from time translation and 
scale-factor symmetry 


The next step is to compute the gauge parameter by 
solving (III. 101) which now reads 


General Covariance essentially guarantees that the cos¬ 
mological evolution equations preserve time translation 
symmetry. For flat space-time, we also have the so 
called scale-factor symmetry. Schematically the differen¬ 
tial equations governing the cosmological evolutions must 
be of the form 

Tl (s3F’ W' SO = 0 with " > 0 “ d ”• 9 - 0 • 

(mi) 

where t is the proper time, which is related to the confor¬ 
mal time via dt = a{t)dr , as usual. In particular, (III.Ill) 
implies that if a{t) = a b {t) is one particular solution of 
the evolution equation, then so is 

a{t) = (1 + a 0 )a b (t - t 0 ) ■ (H.12) 


ao[l + rH\ -toH + ^'+ZH =0 . (11.19) 

The above equation can be integrated to give us 
a 2 (r)£+2 f df a 2 (f)[a 0 (1 + t'H{t))-t 0 H(t)\ + -^- = 0 , 

J T* 

( 11 . 20 ) 

where r* is any fixed time that can be chosen as con¬ 
venient, and we have explicitly included the new dimen¬ 
sionless integration constant bo and introduced M p , the 
reduced Planck mass, for convenience. The integrals can 
be simplified as follows 



dr a 2 (r) (1 + tH(t)) 


1 

2 


a 2 r + 



dr a 2 (r) 


In conformal time the differential equations have a 
slightly different schematic form 


Ti 




= 0. (11.13) 


It is easy to check that for metrics written in the con¬ 
formal gauge this means that the 2-parameter class of 
solutions is given by 


and 



da 1 2 
dr a-— = —a , 
dr 2 


where any further integration constants can be absorbed 
in the bo term. The gauge parameter is thus given by 


£ = -a 0 


ffO) ' 
2 ( r )_ 


bo 

M p a 2 (r) 


+ T 0 ■ 


( 11 . 21 ) 


a(r) = (1 + a 0 )a b [(l + a 0 )(r - r 0 )] 


(11.14) We have defined the function g(r) as 


We emphasize that a b (r) is a specific function describing 
the background which does not contain any free param¬ 
eters (integration constants). 

Let us perturb the general solution around the specific 
background a b (r) to obtain the modes fd s discussed in 
the previous subsection: 


i(r) = (1 + a 0 )a b (T + a 0 r - t 0 ) 

~ a b (r) + aoa b (r) + (a 0 r - T 0 )a(,(r) 

a b {r) 


= a b (j ) 


1 + a 0 + (a 0 r - t 0 ) 


a b (r) 


.(11.15) 


Thus, we have two superHubble modes in the syn¬ 
chronous gauge given by, 


g(r) = f dr a 2 (r) (11.22) 

J T-k 

Finally, one can obtain $ from the gauge parameter 
via (III.71) : 


<F(t) = a 0 (k) 

'Kg 

a 2 

+ b 0 (k ) 

n 

M p a 2 

(11.23) 

or in proper time: 





<F(f) = aoik) 

a 

+ bo(k) 

H ' 

M p a 

(11.24) 


with 


/i(r) = 1 + T—fv = 1 + tT~L (11.16) 

a(r) 

and / 2 (t) = = (H-17) 

a{r) 

associated with the amplitudes (integration constants) 
ao and tq respectively, where we have again dropped the 
subscript “b” for brevity. The corresponding Newtonian 
potential in the synchronous gauge reads 

ip(r) = a 0 [1 + tH\ — ToV. . (11.18) 


g{t) = f dt a(t ) (11.25) 

Jt * 

Note, that we have now included an explicit In¬ 
dependence in ao, bo to emphasize the fact that the am¬ 
plitude associated with the different superHubble modes 
will, in general, vary with the comoving wave number 
k. The ao mode is generally referred to as the constant 
mode while the b o mode is referred to as the growing 
or decaying mode depending on whether the universe is 
contracting or expanding respectively. 
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C. Possible Generalizations 


2. Working with Evolution Equations 


While most of this paper will be devoted to track¬ 
ing the usual superHubble modes associated with the 
ao,h 0 modes, the connection between the superHubble 
evolution and the background solutions discussed above 
is more general and can even be applied to cases where 
the scale-factor/time translation symmetry is not avail¬ 
able, or when there are extra modes present. Here we 
will briefly discuss how our approach and results can be 
extended by looking at certain illustrative examples. 


1. The case of open or closed universe 


If the universe is spatially curved, i.e. open or closed, 
then the scale-factor symmetry is no longer a good sym¬ 
metry of the cosmological equation. For instance, con¬ 
sider the Hubble equation in an open/closed universe in 
the presence of an ideal fluid with energy density, p: 

" 2 = rni( p ~ 5 ) (IL26 > 

Due to the presence of the a dependent curvature term, 
if a(t) is a solution to the above equation, (1 + ao)a(f) is 
no longer a solution. However, there is a way around this 
for an ideal fluid which satisfies the continuity equation 

p + 3TL(l + w)p = 0 , (11.27) 


while, in conformal time, the Hubble equation in the pres¬ 
ence of spatial curvature reads 

" 2 = 3 k { ^- K) (IL28> 

Now, given a solution ab(r),pb(r) of the above cosmo¬ 
logical equations, it is easy to see that 

a(r) = (1 + a 0 )a(T - r 0 ) and p(r) = (1 + a 0 ) _ 2 pb(r - t 0 ) 

(11.29) 

also solves them 3 . Proceeding in a manner similar to 
what was done in the previous subsection, one now ar¬ 
rives at a slightly different result for superHubble modes: 


$(t) 


a 0 (k) 



+ b 0 (k) 


n 

M p a 2 


(11.30) 


We will see later how the above indeed reproduces the 
superHubble fluctuations in a numerical example. One 
may wonder that within GR it should be possible to em¬ 
ploy the same argument for a flat universe as well to 
arrive at (HI.301) instead of (III.231) : is there a discrepancy 
then? The answer obviously is no; in GR for a flat uni¬ 
verse, fig /a 2 , is a constant making (III.301) and (III. 231) 
equivalent to each other. 


3 This trick does not work for more general forms of matter such 
as scalar fields, or if the gravity side of the equation is modified. 


There may be situations where one cannot find solu¬ 
tions associated with the time translation or scale factor 
symmetry or there could be additional modes present in 
the theory. How do we capture the superHubble evolu¬ 
tion in these cases? Well, at least if the cosmological 
differential equations governing the background evolu¬ 
tion is known, then it is possible to derive the perturba¬ 
tion equations that superHubble fluctuations must sat¬ 
isfy. In section IVIIII the appendix, we have illustrated 
this method for GR. Such a prescription may be a lot 
easier than obtaining the perturbation equations by per¬ 
turbing the full gravitational equations. In some case, 
such as in Loop Quantum Gravity, the complete grav¬ 
itational equations may not even be available, but our 
procedure can still yield the equations relevant for su¬ 
perHubble fluctuations. 

Once the perturbation equations are known, one can 
of course solve them, numerically if necessary. What of¬ 
ten is of particular aid in this process is that in many 
phenomenological situations, far back in the past and far 
ahead in the future, the model reduces to GR. The su¬ 
perHubble solutions in this case obviously are known and 
numerical solutions can then be used to connect the su¬ 
perHubble modes of the past with future via “transfer 
matrix”. In a non singular geometry, as long one doesn’t 
use pathological variables it should be possible to ob¬ 
tain these transfer matrices encoding the physics of the 
bounce and turnaround. We will illustrate this procedure 
using a simple bounce-inflation model in section IVll 


III. TESTING THE SUPERHUBBLE 
EVOLUTION FORMULA 

The main aim of this section is to look at GR based 
cosmological scenarios where the perturbations can be 
tracked precisely, at least numerically, and then check 
whether our superHubble formulae, (HL231) and (fTUiOl) . 
holds. In the process, we will also be able to make some 
observations about the subHubble or short wavelength 
fluctuations. 


A. Monotonic Expansions & Contractions 

We are going to first look at simple single fluid mono- 
tonically expanding or contracting FLRW models. Start¬ 
ing with Einstein’s equations 

g ; = t ; , (iii.i) 

one can derive the usual General Relativistic perturba¬ 
tion equations for the Bardeen potential <f>. For a perfect 
fluid (no anisotropic stress) the General Relativistic per- 
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turbation equations for <f> reads 


B. Perturbations Around a Turnaround 


-k 2 $ k -3H(<f>' k + H<f> k ) = ^a 2 Sp k (III.2) 
*l+3H* k + (2H + U 2 )<S> k = ^a 2 Sp k (l II.3) 

The above equations can be combined to yield 

$l+3(l+c 2 s )H$ k + c 2 k 2 <f> k + [2H' + {l + 3c 2 )H 2 ]<$> k = 0. 

(III-4) 

where c s is the sound speed: 



For a single ideal fluid it is related to the equation of 
state parameter, u>: 


1. The toy model 

Tracking fluctuations around turnarounds is easier 
than during bounce primarily because in most cyclic 
models turnarounds occur at energy densities much be¬ 
low the Planck scale where GR is assumed to hold. Thus 
no “new physics” is required to attain a turnaround and 
the evolution equations are well known. Taking inspira¬ 
tion from the cyclic inflation model [52l - l54l [56I . l59j , we 
will here consider a turnaround mediated by a negative 
cosmological constant in the presence of radiation. Ac¬ 
cordingly, the Hubble equation can be written as 

H2 = with p = A(a “ 4 - 1) , (III. 14) 


c 2 = u> where oj = 


P 

P 


(IIII. 41) lias the following super-Hubble limit: 


(HI 6) 


<5>l+3{l+u)U<5>' k + [2'H + (l + 3w)-H 2 ]4> fc = 0. (III.7) 
For a single fluid evolution we further have 

o(r) = \M p r\ q -H = - and W = , (III. 8 ) 

T 


with 


2 

q = - . 

y l + 3w 


(III.9) 


Then the last term in (IIII.4I) drops out and we are left 
with 




6 ( 1+ ^ 

r( 1 + 3w) 


= 0 


(III.10) 


In the above expression, r > 0 represents an expanding 
phase while r < 0 corresponds to a contracting universe. 
(IIII. 101 ) has two well-known solutions, the constant mode 
and the decaying/growing mode: 


$ = C + 


D 

\M p t\ 2c i +1 


(III.ll) 


Let us now verify that (III.231) indeed reproduces the 
above two modes. Let’s specialize to the expanding uni¬ 
verse, the calculations for the contracting universe follow 
similarly. For (IIII. 81) the function g is given by 


9 = j dra 2 {T) = M p 1 ( M p r) 2q+1 (III. 12) 

(III.231) then yields 


<F(r) = 


3ao(l T cc) 2hq 

5 + 3w + (1 + 3uj)(M p t) 2c > +1 ’ 


(III.13) 


which is indeed of the form (IIII. 111) . 


We have chosen our conventions such that the 
turnaround occurs at a = 1. A is a free parameter and 
governs the time scale of the turnaround. For numerical 
simulations it becomes imperative that we work with the 
Ray Chaudhury (RC) equation 4 : 


a 

a 


1 

6 M 2 


{P + 3 p) 


where 


P = Pr + A = A 



(III.15) 


(III.16) 


Also, in order to compare our numerical results with the 
analytical expressions of (III.231) we have to work with 
conformal time. The Hubble and RC equations read 


= “ d T = 6SJ ( '’- 3f ’>- < IIU7 > 

Let us next look at the perturbation equation, (IIII. 41) . 
Since the cosmological constant does not contribute to 
any fluctuations, the fluctuations in pressure and energy 
density only comes from radiation, and therefore we have 


(IIII. 41) then can 
and RC (IHLITI) 

+ 4 H$' k 


2 _ 4 
° s ~ 3 


(III.18) 


be simplified using the Hubble (IIII. 141) 
equations to yield 


1 .9 t 4A 

- -( 

3 3M 2 


= 0. 


(III.19) 


4 Numerically, the Hubble equation can only be solved separately 
for the expanding and contracting phases as one needs to take 
its square root before it can be implemented in a code. Thus it 
is not possible to track the evolution of the background or the 
perturbation through bounces and turnarounds using the Hubble 
equation. The Hubble equation, however, acts as a constraint 
determining the initial condition for a '. 
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2. Numerical Analysis 


It is rather useful to work in units of the time-scale 
associated with the bounce: 

[W 

T=\j -jf (III.20) 

Accordingly, the relevant dimensionless variables are 

t= 1 ; H = m= a,nd k = kT (III.21) 
/ a dr 

The evolution equation now becomes 

K + 4 UA>' k + - ^a 2 $ fc = 0 , (III.22) 

where the prime now denotes differentiation with respect 
to the t variable. The above equation can be solved 
in conjunction with the acceleration equation (IIII.17I) . 
Note, that the initial condition for a' can be determined 
from the rescaled Hubble equation as follows: 


da I a 4 T 2 p 11 — a 4 

df = + \J 3M p 2 = V 3 


(III.23) 


comparative “yard-stick” in determining what consti¬ 
tutes a superHubble fluctuation. For instance, exactly 
at the turnaround point, H = 0, and therefore techni¬ 
cally, all modes are subHubblc! Intuitively however, we 
know that such a criteria is not physically sensible; for 
instance, for sufficiently long wavelengths the time their 
wavelengths remain smaller than H is too short to even 
establish a causal connection over the wavelength itself. 
Thus the “subHubble” phase should not have any effect 
on evolution for these modes, and the modes should es¬ 
sentially follow the superHubble trajectory. This is ex¬ 
actly what our numerical simulations indicate: Modes 
with sufficiently small k -C T -1 follows the analytical su¬ 
perHubble trajectory (III.231) . In other words, the correct 
time scale that one should compare the physical wave¬ 
lengths of fluctuations to determine whether the mode re¬ 
mains “superHubble” is the turnaround time scale, given 
in this model by T. We will see that the similar argu¬ 
ments also hold for bouncing cosmologies. This is an 
important result in the context of cyclic/bouncing cos¬ 
mologies as it greatly simplifies tracking long wavelength 
fluctuations. 



FIG. 1: Numerical evolution of <f> for kT = 0.2, 0.4, 0.6, 0.8, 
1.0, 1.2, 1.4 (red, orange, ochre, green, blue, indigo, purple) 
during a turnaround. The black curve corresponds to the 
analytical result for the superHubble evolution. Evidently, 
the colored curves approach the black as k —> 0. The initial 
conditions for the above graph corresponds to o(0) = ao = 0.1 
and 4>'(0) = 0. 

Fig. [I] depicts the evolutions of <f> for different values 
of the comoving wave number k = kT. The red curve 
denotes the analytical result (III.231) and as is evident, 
the blue curves approach it in the k —> 0 limit. Firstly, 
this corroborates the validity of our separate universe 
approach to track the perturbations during transitions 
between expansion and contraction, and in particular the 
formula (IIL23I) . 

Secondly, it supports the view point [53 l. [60| that the 
Hubble parameter is not always the most appropriate 


3. SubHubble Fluctuations 

What can we say about the evolution of the subHub¬ 
ble fluctuations? Unlike the superHubble fluctuations 
where the spatial fluctuations did not matter because 
their wavelengths were too large as compared to the cos¬ 
mological evolution, the same cannot obviously be argued 
for subHubble fluctuations. Moreover, often QG theories 
involve higher derivatives, which would necessitate keep¬ 
ing track of higher powers of k in the evolution equations. 
The situation is complicated and model dependent, and 
to our knowledge no general approach to tracking sub¬ 
Hubble fluctuations have been discussed in the previous 
literature. In Fig. [2] we have plotted the evolution of 
4> for different subHubble wavelengths and they oscillate 
with time-varying amplitudes. Can we find any general 
trend? 

It is known that for monotonically expanding or con¬ 
tracting GR models with a single fluid or a scalar field 
although, the amplitude of <f> oscillations fluctuate, the 
oscillations in the so-called “canonical” variable, y, have 
constant amplitude f58). This is because the subHub¬ 
ble wavelengths are so short as compared to the cos¬ 
mological time scale that one can find a canonical vari¬ 
able where the cosmological dynamics almost freezes and 
the corresponding action resembles that of a free-held 
in Minkowski space-time. What happens to y during a 
turnaround? 

For a cosmology dominated by ideal fluids the canoni- 
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FIG. 2: Top: Numerical evolution of $ for subHubble modes: 
kT = 1, 10, 25, 100 (red, green, blue, purple) during a 
turnaround. The black curve corresponds to the analytical 
result for the superHubble evolution. Bottom: Numerical 
evolution of x f° r the same modes. Evidently, the amplitude 
of fluctuations in x stays a constant for the subHubble modes 
except at the bounce point where it is ill-defined. The initial 
conditions for the above graphs correspond to a(0) =0.1 and 
<F'(0) = 0. 

cal variable is defined via [H) 

X = (III.24) 

wl,h « = t + 3(T^(l + 4 ) -< IIL25 > 

and 2 = a\j 1 — ■ (III.26) 

Clearly, \ is not a good variable exactly at the 
turnaround point since z diverges at this point. The 
same happens during a bounce as we will see in sec¬ 
tion ED Sometimes the divergence of \ is taken as ev¬ 
idence that a bounce is inherently unstable, but as seen 
in Fig. [T] and El the metric is completely well-defined and 
the perturbations can remain well under control. In sec¬ 
tion ED we will see the same phenomenon repeating dur¬ 
ing a bounce. This is a clear demonstration that bounces 
and turnarounds are not inherently unstable and the cos¬ 
mology can remain true to the FLRW assumption of a 
homogeneous universe. 

What is perhaps more interesting is the observation 
that the amplitude of subHubble fluctuations in \ does 
not change during the course of the bounce, see Fig. [2] 


(bottom). Thus our intuition based on monotonic cos¬ 
mologies in GR about x oscillations seems to continue to 
hold. If such a result could be proven for more general 
situations, it would be a tremendous step in our under¬ 
standing of perturbative evolutions, and in section IVI1 
we will argue why there is good reason to be optimistic 
about this trend. 


C. Perturbations Around Bounce 


1. Toy model with gas of branes 


In this subsection, the plan is to study how perturba¬ 
tions evolve around a bounce. For the purpose of illustra¬ 
tion, we are again going to consider a toy model that can 
produce a bounce within the confines of the General Rel¬ 
ativity. We are going to look at a closed universe model 
containing a perfect fluid satisfying —1 < w < —1/3. 
In this case the background cosmology gives rise to a 
bounce inflation scenario where a fast contraction is fol¬ 
lowed by an accelerated expansion mediated via a non¬ 
singular bounce. Inflation is of a power-law type which 
approaches the exponential case in the limit w —> — 1. A 
particular example of an ideal fluid with such negative 
equation of state parameter is a gas of co-dimension 2 
branes (or domain walls), which has an equation of state 
to = —2/3. It is also possible to mimic such equation of 
state parameters with a gas of strings appropriately cou¬ 
pled to scalars such as the dilaton, see for instance [6li |. 
We are however not really interested in the phenomenol¬ 
ogy of such models; our purpose is to test our formalism 
for tracking superHubble fluctuations in bouncing back¬ 
grounds where explicit numerical calculations can be per¬ 
formed. 

With this in mind, let us proceed to obtain the differ¬ 
ential equations that govern the evolution of the back¬ 
ground and the perturbations. The total “effective” en¬ 
ergy density and pressure that enters the Hubble equa¬ 
tion (MI. 141) and the RC equation (IIII. 171) can be written 
as 

P = Pf~ P -\ = P.(a- 3(1+W) - a“ 2 ) (III.27) 

ar 

P = Pf + ~ = P. (w«- 3(1+u) + ^a- 2 ) C1II.28) 


where we have chosen a convention such that the bounce 
occurs at a = 1. 

The perturbation equation for the Bardeen potential 
needs to be modified to account for the curvature of the 
universe and according to our conventions reads [Hsj ] 


<■+3(1 + c 2 s )H<S>' k +c 2 s k 2 <t> k 


2 H' + (1 + 3c 2 ) 




0 .(III.29) 


The fact that it is negative simply indicates the instabil¬ 
ity of the system at short distances, but here we are only 
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interested in the long wavelength fluctuations. In Fig. 0 
we have illustrated the behavior of $ for different values 
of k for the domain wall case corresponding tow = —2/3. 

In Fig. [3] we have plotted versus the re-scaled con¬ 
formal time r — > t/t,, where r. = M|>/ p, is the bounce 
time scale. Again, we see that as k —> 0, the evolution ap¬ 
proaches the analytically derived superHubble evolution 

mm- 
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FIG. 3: Numerical evolution of <f> for fcr. = 0.05, 0.1, 0.2, 
0.3, 0.5, 0.8, 1.0 (red, orange, ochre, green, blue, indigo, pur¬ 
ple) during the bounce. The black curve corresponds to the 
analytical result for the superHubble evolution. Evidently, 
the colored curves approach the black as k — > 0. The initial 
conditions for the above graph corresponds to a(0) = 10 and 
$'( 0 ) = 0 . 


IV. APPLICATION I: TRANSFER MATRIX 
AROUND BOUNCE 

In the previous sections we have derived analytical ex¬ 
pressions for the time evolution of the two superHubble 
modes, and verified their validity in cosmological models 
involving monotonic expansions or contractions as well 
as when there is a transition between the two phases. 
We performed these checks in GR based models where 
an exact numerical analysis could be performed. In the 
next couple of chapters our goal is to apply our derived 
superHubble formula (III.231) to understand the evolution 
of perturbations for certain classes of nonsingular cos¬ 
mological models. On one hand, this will serve as an 
illustration of how one can apply our formalism and re¬ 
sults to different cosmological models, and on the other, 
it will shed insights into these specific scenarios that we 
hope will be useful in building viable and consistent cos¬ 
mologies. 

In cosmology, we are most often interested in calculat¬ 
ing the “primordial” spectrum of fluctuations that seeds 
the temperature anisotropies observed in the CMB. The 
usual algorithm is the following. One starts with a choice 
about the initial conditions (IC’s) for the fluctuations in 
the subHubble phase. Traditionally, one uses the Bunch- 
Davis vacuum to obtain the IC’s. Although other vacuum 


choices [62] as well as statistical thermal initial conditions 
have been explored [22|, [63|, HH, fortunately, the under¬ 
lying algorithm for calculating the primordial spectrum 
does not really depend on the IC prescription. Typi¬ 
cally, the fluctuations are tracked as they pass from the 
subHubble to the superHubble phase at which point the 
perturbations are “expected to freeze”. In this paper, 
primarily we are concerned with the superHubble evo¬ 
lution that ensues. Do perturbations really freeze, even 
when there are cosmological transitions such as bounces 
and turnarounds? How is the “initial” spectrum that is 
generated while Hubble crossing (possibly during a con¬ 
traction phase) related to the “primordial spectrum” in 
the expanding branch that is suppose to seed our CMB 
fluctuations? 

Of course, we already have an analytical formula for 
the superHubble fluctuations, but to address these ques¬ 
tions it is often more convenient to obtain a “transfer 
matrix” that relates how the amplitudes of fluctuation of 
the decaying and constant mode of the expansion phase 
is related to the growing and constant mode of the con¬ 
traction phase when the expansion and contractions are 
connected via a nonsingular bounce or turnaround. In 
this section, we are going to illustrate how our formalism 
can be used to answer these questions within the frame¬ 
work of a simple bouncing cosmology. 


A. A Simple Bounce Cosmology 

Let us consider a simple but consistent example of a 
bouncing universe dominated by a single fluid with w > 0, 
so that, unlike the brane-gas scenario discussed in sec¬ 
tion IIII Cl the underlying thermodynamic system is sta¬ 
ble. The price we have to pay is that we can no longer 
stay within the framework of GR as various no-go the¬ 
orems exists thatprecludes bounce in GR with regular 
(w > 0) matter [J-Q. Thus, we envision having a con¬ 
tracting GR phase that lasts approximately till Planck- 
ian energy density is reached. This is followed by a QG 
bounce with a Planckian cosmological time scale that 
eventually leads to a GR expanding branch once the en¬ 
ergy density has diluted to sub-Planckian values. Now, 
during such an evolution, there are going to be some 
modes which are so short that even at the bounce point 
their wavelengths are shorter than the cosmological time 
scale, Mp 1 ~ 10“ 42 s. These modes remain subHubble 
throughout the course of the evolution and we will not 
discuss their evolution any further. All other modes exit 
the Hubble radius some time in the contracting GR phase 
and remain superHubble during the course of the bounce. 
Let us remind ourselves that since the time scale of the 
bounce is ~ M” 1 , all the modes whose physical wave¬ 
length is larger than the Planck length near the bounce 
behave as super-Hubble modes. 

Now, as discussed in section lHI Al in the GR regime the 
super-Hubble fluctuations are known to be of the form 
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FIG. 4: Schemtic Hubble diagram for an u> > 0 universe. 
Black curve represents the cosmological time scale, while the 
colored curves correspond to the physical wavelengths for spe¬ 
cific comoving modes. All the modes start out sub Hubble in 
the GR phase. Most of them exit the Hubble radius during 
the contracting phase, and re-enter the Hubble radius in the 
expanding branch. Some very short modes (purple curve) can 
always remain sub Hubble. 


tnnn 


= c± ± ■ (IV ' 30) 

The — and + indices referring to the super-Hubble phases 
in the “past” and the “future” branch respectively. To 
gain insights into the evolution we want to find the trans¬ 
fer matrix relating C+, D + with C_, D-. For instance, if 
the matrix is found to be diagonal then we can conclude 
that there is no mixing between the modes; the growing 
mode during contraction goes over to the decaying mode 
during expansion, while the constant modes on the two 
sides of the bounce map to each other. Another impor¬ 
tant issue has to do with whether there is any amplifi¬ 
cation/reduction in the amplitude of fluctuations as we 
go through the bounce. As a final point, we should also 
check whether the fluctuations can remain small during 
the entire bounce phase for us to justify using perturba¬ 
tion theory. 

With these goals in mind, let us begin by choosing a 
coordinate system such that r = 0 corresponds to the 
bounce. Then, in order to obtain the transfer matrix, 
all we have to do is to expand the solution (III.231) that 
is valid in the entire super Hubble phase as a series in 
(M p r) _1 for large times 5 . The leading and next to lead¬ 
ing order terms should then correspond to the two super 
Hubble modes in the GR phase. This way one can de¬ 
termine C±,D± in terms of bo,ao , and subsequently the 
transfer matrix. 

Now, we know that as long as |r| M~ x we can trust 


5 This is assuming that non-GR effects appear at Planck scale 
which is then also identified as the bounce scale. If some lower 
mass scale effects (such as may be from string theory) are present 
then one has to look at appropriately larger times. This however 
doesn’t change the basic algorithm proposed here. 


the GR evolution and accordingly 

a E {j) « e±\M p r\ q , 


(IV.31) 


where the +, — indices again indicated the regions, 
M p t 1 and M p t 1 respectively. For a symmetric 
bounce e+ = e_, but not otherwise. 

Next, let us define 

g = f a 2 (f)df (IV.32) 

Jo 

While evaluating the integral for positive r, we can break 
up the above integral into two parts: [0,r.] and [r.,oo), 
with r. M~ x . In this case g(r) can be evaluated as 
follows: For r > r. 


g(r) = dr a 2 (r) + dr cr (r) 


dr a 2 (r) + dr a 2 E (r) 


dr [o 2 (t) - a|(r)] + dr a 2 G (r) 


= 9+ + [ dr a|(r) 
JO 


(IV.33) 


The second term can be integrated straightforwardly to 
give us 


g(r) =g+ + e+ Wp T ) 2q+ * for > (IV.34) 

’ y+ M p (2q + 1) ~ V ; 

We note that since for large r, a(r) —> a^(r), and g + 
tends to a constant as one increases r., as it must for the 
approximation to be consistent. For r < —r, one can 
perform a similar analysis by breaking the integral as 



leading to 


g( T ) 


g- 


e_(—A/ p r) 2g+1 ' 
AI p (2q +1) 


for r < —r. , (IV.35) 


where 


g± = ± I dr [u 2 (t) - a|(r)] (IV.36) 

Jo 

Again, for a symmetric scale factor we should have g + = 
g-. To summarize, quite generically we have 


a(r) 

-A e±\M p T\ q 

(IV.37) 

n 

£ 

T 

(IV.38) 

g 

[ e ± |M p r| 2 ^l 

r + m p ( 2,+i)J • 

(IV.39) 
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We are now in a position to evaluate $ for large times: 


$ ± (r) = a 0 


2q±l 


- 1 


q(a 0 M p g + ± bp) 
\M p t\ 2 i +1 


. (IV.40) 


We now explicitly see that for |r| > r # , the first and 
the second terms in (IIV.40I) correspond to the constant 
mode and the growing/decaying mode respectively. By 
matching the above asymptotic expansion with the GR 
super-Hubble solutions (1IV.30I) . we can find the transfer 
matrix that we are after. We have 


C± — ao 


e±9 


- 1 


_ 2q + 1 
D± = q(a 0 M p g± ± b 0 ) 


(IV.41) 
(IV.42) 


One can rewrite the above equations in a transfer matrix 
language: 


D \ 



/ q(e+~ 2 ) —1 

= ( q( e - —2)—i 

l <j(2g+l)M p (g + +g_) 

\ g(e_-2)-l 

(IV.43) 

Several comments are now in order: Firstly, it is now 
clear that while the modes do mix, the constant mode 
in the expanding branch does not get any contribution 
from the growing mode in the contraction phase. This 
is a crucial result for the purpose of model building in 
bouncing and cyclic models. For instance, it corrobo¬ 
rates the result (57j, that in a single fluid scenario, at 
least, the scale-invariant perturbations of the growing 
mode in ekpyrotic cosmology cannot be transferred to 
the constant mode in the expanding branch that is rel¬ 
evant for CMB. Secondly, we see that for a symmetric 
bounce C+ = C-, i.e. , there is no amplification in the 
constant mode. However, if the bounce is not symmet¬ 
ric, i.e. e_|_ ^ e_, the amplitude may change. In a single 
fluid universe it is hard to imagine how one can realize 
an asymmetric bounce, but our results strongly suggest 
that amplification may be possible in a multi-component 
asymmetric bounce. 

Finally, let us discuss the general behavior of the super- 
Hubble perturbations near the bounce. Let us look at the 
bo mode first. We see that in the GR contracting phase 
the mode grows very fast ~ \M p T \~ (y2q+l ' 1 . However, as 
the universe enters the QG phase around r ~ M” 1 , 
T-L —> 0, and the scale factor reaches a minimum, so that 
approximately one expects its largest value to be given 
by 'H/MpO? ~ 0(1). For the validity of the perturbative 
analysis then 6o 1, but importantly, the perturbations 
do not grow unboundedly during the bounce. 

What about the ao mode? Clearly it begins as a con¬ 
stant in the infinite past. According to our conventions, 
while \H\/a 2 keeps growing in the contracting phase \g\ 
decreases, so that the overall growth of the ao mode is 
slightly weaker than the bo mode. Exactly at the bounce, 
again g = H = 0 while a is at its minimum (but finite) 
value. Thus again, the perturbations remain bounded 
and attains its maximum value of \gT~L/a 2 \ ~ 0(1) around 


r ~ M p . To summarize, as long as ao, bo 1, although 
both the superHubble modes grow towards the bounce, 
they can remain small enough for perturbation theory to 
remain valid, and the bounce scenario to be consistent! 


B. Illustrative Examples 


Our formalism is particularly useful in situations where 
the evolution equations for the perturbations are either 
not known or too complicated. Let us therefore illustrate 
our algorithm in one such example; bouncing cosmology 
arising in Loop Quantum Cosmology where although the 
background equations have been studied in great depth, a 
complete understanding of the perturbations is still lack¬ 
ing. We will start with the modified Hubble equation 
that appears in LQC [29j[3C|: 

(IV44) 

where p . ~ M p . Thus at low energy densities, p <C M p , 
we recover the usual Hubble equation. 

For an ideal gas, this gives us 

D fl 

H 2 = * M2 (a n ~ 1) with n = 3(1 + u>) , 
which can be recast into 


da n 

dt 


= n. 


I p,(a n — 1) 
3 M 2 


A straight forward integration yields 


a — 1 + 


n 2 t 2 




1/n 


12 M p J 


(IV.45) 


(IV.46) 


For large values of t we clearly recover the GR limit. 

We are now going to expand the LQC perturbed so¬ 
lution in the small variable x = l/\M p t\ to obtain the 
superHubble limit in the GR phase. Specifically, for the 
background solution (IIV.46I) . one finds 


a(t) 

H 


zz2l ( , x 

OL n X n I 1 H- 

a 

± 2 aMpX 


2 \ n 


a n x 


i+-) 

na J 


9 = 


n(l + x 2 /a ) n 

i i : 

2 ,_ n’ 2’ 


2 

—aM„x 


x 2 

1 - — 

a 


± 1(1 13 a 

2 Ai ( — 


(IV.47) 


where 


_ n 2 p. 

~ 12M£ ’ 


(IV.48) 
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and 2-Pi is a hypergeometric function with the limiting 
expression 


2F1 


/1 13 a\ 


= x 




2r(-i) 


+ o(^) 


+ 


m 


n + 2 


C(a 


(IV.49) 


as x —>■ 0 or M p —» 00. Thus the dominant terms in the 
large t limit are given by 


9 = 


±1 



— n + 2 1 . 
+ x n a n 



(IV.50) 

To compare our results with previous analysis we will 
need to know how the conformal time is related to phys¬ 
ical time t in this model. We choose 


t = 



(IV.51) 


which for the GR limit gives us 

\M p t\ = [(1 + q)\M p t\]T*i (IV.52) 

=► {l + q)^\M p t\^ = \M p r\ 1+2q (IV.53) 


It is easy to check that a, g indeed reduces precisely to 
the generic form enumerated in (IIV.37I) and (IIV.39I) with 


9+ = 9- 
and e + = e_ 


2M P T (-1) 
a™ (1 + q) " . 


(IV.54) 
(IV.55) 


Accordingly, the transfer matrix in an LQC bounce is 
given by 


T = I g(2g+i)y/fr(-i-£) 

r (-^)[9(«"( 1 +9)“"-2)-l] 


0 

-1 


(IV.56) 


V. APPLICATION II: STABILITY OF THE 
EMERGENT CYCLIC INFLATION MODEL 

A. Emergent Cyclic Inflation: a brief Review 

As an application of the superHubble tracking formula 
dirm in this section we will look at the evolution of 
superHubble fluctuations in the emergent cyclic infla¬ 
tion (ECI) scenario, a new inflationary paradi gm that 
has been proposed in a series of articles @, f5lt455l @. 
In the presence of a bounce mechanism, ECI hopes to 
realize a geodesically complete non-singular inflationary 
space-time via a series of asymmetric short cycles. Our 
primary aim is to check whether such a cosmology can be 
stable. For an overview on the model and its motivations 
see fHHj], but here is the ECI paradigm in brief: 


• Cyclic Emergent Phase: The universe “begins” in 
a quasi-periodic phase of oscillation. By quasi- 
periodic we mean that in each cycle the universe 
grows a little more than it contracts. These oscil¬ 
lations become more and more symmetric and pe¬ 
riodic as we go back towards the past infinity [6|. 
The cosmology is very similar to the “emergent uni¬ 
verse” scenario advocated in ill- The emergent 
phase, as envisioned in @, works if the universe 
is closed and contains different species of matter 
which are not always in thermal equilibrium with 
each other so that in each cycle a little bit of en¬ 
tropy is produced giving rise to the slight asymme¬ 
try in the cycles. While this asymmetric growth is 
required for the universe to eventually be able to 
exit the cyclic emergent phase, if this asymmetry 
persisted all the way to past infinity, that would 
lead to a geodesically incomplete universe with the 
amplitude of scale factor fluctuations vanishing as 
t —> — 00. However, it was first argued in d and 
later numerically verified in [55 ] that, in the cyclic 
emergent model as one goes back in the past the 
entropy production becomes less and less in every 
cycle and the universe asymptotes to a periodic os¬ 
cillatory phase as t —» —00. Such a space-time 
geometry is obviously geodesically complete. 

In the cyclic emergent phase, the turnarounds oc¬ 
cur when the negative curvature density associated 
with the closed universe cancels the total matter 
density, while for the bounces to occur one needs a 
separate mechanism in place. 

• Cyclic Inflation Phase: In the ECI scenario pre¬ 
sented in [53I, [HU [56}, as the universe grows slowly 
in the cyclic emergent phase, there comes a mo¬ 
ment when the negative potential energy 6 (of vari¬ 
ous scalar fields that may be present) becomes more 
important than the negative curvature density. At 
this point the universe transitions to a cyclic infla¬ 
tion phase where the turnarounds now occur when 
the matter density is canceled by the approximately 
constant negative potential energy. One can easily 
check that now in every cycle has the same time 
period, and the volume of the universe grows by 
the same factor. The latter occurs because the in¬ 
teraction between different particle species tend to 
naturally increase the entropy of the universe by 


String theory, the leading candidate for a consistent theory of 
quantum gravity, naturally predicts the existence of negative 
energy vacua, and it has been quite a challenge to find ways 
that may lead to positive vacuum energies in the string the¬ 
ory/supergravity framework [65l , |66) . In general, if one looks 
at the potential energy coming from all the moduli in any fun¬ 
damental theory, one would expect to have both negative and 
positive potential regions, possibly with several local minima’s 
dispersed liberally. In String theory this picture is often referred 
to as the “landscape” [67) . 
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the same factor: 

= 1 + 3k , (V.57) 

where S n denotes the entropy of the “nth” cycle. 
Since the total entropy is proportional to the vol¬ 
ume, this means that for small k. the scale-factor 
(at say the bounce point) in consecutive cycles in¬ 
creases by a factor (1 + k). Over many cycles the 
space-time resembles the inflationary growth [52j |! 

• Graceful Exit: In [54] it was shown that the Cyclic 
inflationary phase can end if scalar field dynamics 
can classically propel the universe from a negative 
to a positive potential energy region. While this 
depends on several factors, such as how steep is 
the transition potential from the negative to posi¬ 
tive potential regions, and what is the phase of the 
universe when the transition begins to occur, etc., 
it was found in [54[ that for most “natural” model 
parameters, it is much more likely for the universe 
to make the exit than not. 


the universe has to undergo an infinite set of cycles be¬ 
fore reaching our current phase of monotonic expansion, 
and therefore encounters an infinite number of Planckian 
density bounces. Moreover, in the Cl phase, after the co¬ 
moving perturbations exits into the superHubble phase, 
the universe still performs several oscillations 7 . Do the 
spectrum and amplitude of the perturbations relevant for 
CMB change during this process? Let us try to answer 
some of these questions with the help of the formalism 
we have developed in this paper. 

We will start with the emergent phase where the uni¬ 
verse is undergoing approximately a periodic cyclic evo¬ 
lution. A simple model that provides such a periodic evo¬ 
lution is provided by a LQC mediated bounce in the pres¬ 
ence of a negative cosmological constant and radiation. 
This is essentially an extension of the scenario discussed 
in section UlI Bl The LQC correction (IIV.44I) provides a 
bounce when p = p., while the turnaround occurs when 
p r = A. This gives us a periodic cyclic evolution of the 
universe much like one envisions in any emergent cosmol¬ 
ogy 8 . For numerical purposes, one again has to resort 
to the acceleration equation which reads 


• “Our” Expanding Universe: After the graceful exit 
the potential energy becomes positive and there¬ 
fore, after one last bounce, the universe cannot 
turnaround any more. It has no other choice but 
to continue to expand monotonically which is where 
presumably we now find ourselves in. No reheat¬ 
ing is required in the ECI paradigm as radiation 
remains the dominant energy density component 
throughout the entire evolution. 

To summarize, the ECI paradigm provides a nonsin¬ 
gular geodesically complete inflationary cosmology, does 
not require any reheating (which has it’s own chal¬ 
lenges [68[), and interestingly predicts distinctive sig¬ 
natures in the form of small oscillatory wiggles in the 
power spectrum (53j. With the possible exception of the 
monopole problem the ECI can address all the old cos¬ 
mological puzzles as in standard inflation. Last but not 
the least, the Cl mechanism provides new avenues of con¬ 
structing phenomenologically viable inflationary models 
where the “initial” potential energy is negative. Keeping 
this in mind, here we want to address some of the ques¬ 
tions that is pertinent to ascertain the viability of the 
ECI models. 


B. The Homogeneity Problem in Emergent Phase 

A key problem that often plagues cyclic and bounc¬ 
ing models is whether the assumption of an approxi¬ 
mately homogeneous and isotropic cosmological metric 
holds during the bounce? The concern is well justified 
as it is known that both anisotropy and inhomogeneity 
increases as the universe contracts to high densities. The 
problem is more profound in a model like ECI where 


a 

a 


1 

Ml 


g(p + 3p) 


— (3p + 2p) 

P» 


(V.58) 


where p and p refers to the total energy density and pres¬ 
sure given by (IIII.14I) and (IIII.16D . 

(IV.58I) can be solved straightforwardly, and then the 
two superHubble mode functions associated with ao,bo 
can be plotted, see Fig. [5] (top). The results are rather 
interesting: we observe that while the bo mode oscillates 
with constant amplitude, the ag mode starts to grow. 
This can be understood as follows: The expression for 
the b o model contains 7i and a, both of which are peri¬ 
odic functions and therefore it is no surprise that the bo 
mode remains periodic. In other words, if we start in a 
“patch” of the universe where bo amplitude is small, it 
remains small during the emergent phase. The ao mode 
however, contains the function g which is an integrand 
over a positive definite function and therefore can only 
increase. This results in the growth of the superHubble 
fluctuations in consonance with some earlier findings [25l — 

H3- 


7 The number of oscillations depend on the parameters of the 
model but for certain favored parameter sets, it was found that 
the universe undergoes around 0(1O 3 ) cycles after the modes 
that we see in the CMB enters the superHubble phase [531 . 

8 In the emergent-cyclic model the cycles are actually not exactly 
symmetric. Due to entropy production, in each cycle the universe 
grows a little bit more than it contracts. This is important for 
the universe to be able to exit the cyclic-emergent phase to the 
cyclic inflation phase. This slight asymmetry however, should 
not affect the evolution of the perturbations in any significant 
way and therefore we have decided to exclude the entropy gener¬ 
ating mechanism in our analysis. For more details on this please 
see [55)| . 
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This is a rather important and interesting result for 
the ECI model as it suggests that only patches which are 
smaller than the cosmological time scale of the bounce 
can remain smooth and homogeneous; fluctuations with 
wavelengths longer than this time scale becomes large. 
For a Planckian bounce, this means starting out with 
Planck size patches. 

Addressing whether the amplitude of the subHubble 
(subPlanckian, in our case) fluctuations can remain small 
or not, is beyond the scope of our paper, but in a radia¬ 
tion dominated universe radiation pressure is expected to 
be able to prevent the growth of small scale fluctuations. 
In other words, Jean’s instabilities are only expected to 
grow at wavelengths that are larger than the cosmological 
time scale. This is what we will also see in the following 
section in a scalar field driven bounce inflation model. 
Thus, a plausible picture that emerges in the ECI sce¬ 
nario is that of a “patchy” universe, where some of these 
Planckian patches can be homogeneous and isotropic 9 
that may eventually lead to a universe such as ours after 
a cyclic inflationary phase. 




C. SuperHubble Fluctuations in Cyclic Inflation 
Phase 


We have just witnessed that as long as the cyclic evolu¬ 
tion is periodic the superHubble modes grow and there¬ 
fore one can only envision a smooth isotropic patch of 
the size of the bounce time scale. We now want to ask 
what happens to the superHubble perturbations if the 
universe starts to grow as envisioned in cyclic inflation. 
The detailed dynamics of the asymmetric cyclic growth 
depends on the specific model of cyclic inflation, but the 
advantage of the formalism that we have developed is 
that we only need to know the background cosmological 
evolution in order to track the long wavelength fluctu¬ 
ations. For the purpose of illutration, we will therefore 
look at a simple cyclic inflation background: 

a(t) = e xt [l + esin(uA)] with e <C 1 and uj A . (V.59) 


In Fig. [5] (bottom) we have plotted the two mode func¬ 
tions for the above cyclic inflation background. We ob¬ 
serve that one of the mode decays while the other mode 
oscillates with a constant amplitude. This can be under¬ 
stood easily by inspecting (III.241) : In the Cl phase, H is 
periodic but a grows. This means that the bo mode now 
decays as ~ 1/a. The a o mode, however, contains an 
additional g function whose growth oc a(t). This cancels 
the 1/a decay to lead to a constant amplitude! 

This is a very important result: Firstly, this means that 
the superHubble fluctuations do not grow and hence once 
we start with a homogeneous and isotropic (Planckian) 
patch, the superHubble fluctuations won’t spoil the un¬ 
derlying FLRW behavior. Secondly, the ao mode which 
ultimately becomes the constant mode in the final ex¬ 
panding background, seems to be able to retain its pri¬ 
mordial amplitude of fluctuation through the course of 
the cyclic evolution. This is what was assumed in [53] 
in order to successfully explain the nearly scale-invariant 
spectrum in CMB, and we have now verified the conjec¬ 
ture. To summarize, our analysis gives credence to the 
ECI model with a twist that one can only assume a ho¬ 
mogeneous and isotropic patch of the size of the Planck 
length (or bounce scale) to begin with, but that is a per¬ 
fectly acceptable scenario. 


FIG. 5: The evolution of the ao mode (red) and the bo mode 
(blue) in the emergent cyclic (top) and the cyclic inflationary 
(bottom) phase. 


9 The anisotropies are expected to grow as ~ 1/a 6 . Since in the 
emergent cyclic phase the cycles essentially become periodic at 
past infinity, the scale factor at the bounce point asymptotes to 
a finite value. Unlike inhomogeneities, the anisotropies therefore 
don’t grow and can remain under control. 


VI. APPLICATION III: BOUNCE INFLATION 
A. Background Cosmology 

Our final study involves a simple realization of the 
“Bounce Inflation” scenario where the usual inflationary 
cosmology is preceded by a nonsingular bounce. The 
bounce provides past geodesic completeness to the infla¬ 
tionary space-time and envisions a cosmological evolution 
bereft of the Big Bang singularity and conceptual issues 
related to the beginning of time. Phenomenologically, it 
was pointed out recently fl5l ] that the bouncing phase 
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may be able to account for the anomalous blue tilt that 
has been observed in the low-Z part of the CMB spec¬ 
trum num. Thus a study of perturbations in such a 
model is interesting in it’s own right, but our primary 
interest lies in checking how well the formalism we have 
discussed in tracking superHubble perturbations can be 
applied to this scenario and whether we can develop any 
intuition when it comes to the subHubble fluctuations. 

To this end let us consider a standard scalar field driven 
inflationary paradigm evolving according to the usual 
Klein-Gordon equation: 

ij> + 3= 0 , (VI.60) 

acp 

which in conformal time becomes 

0" + m<\>' + = 0 . (VI.61) 

acp 

The energy density and pressure of the scalar held are 
given by 

P0 = + V(</>) = \(^) + V{4>) , (VI.62) 

IV = - V{<t>) = \ (£) - Vtf) . (VI.63) 

For the purpose of illustration, we are going to assume a 
simple slow-rolling potential, 

V(<f>)=M 4 (v o + tqAj, (VI. 64) 

with ui < 1, which can drive a near exponential inflation. 
For a spatially hat universe the space-time, however, is 
geodesically incomplete. This problem can be avoided, if 
instead, we have a closed universe. In this case, the total 
“effective” energy density and pressure that enters the 
Hubble equation (IIII.14D and the RC equation (IIII.17I) 
are given by 

M 4 

P = P<t> - — (VI.65) 

or 

1M 4 , 

V = IV +o— (VI.66) 

O CL 

The universe starts in a phase of contraction where for 
a sufficiently hat potential the scalar energy remains al¬ 
most a constant while the negative curvature density in¬ 
creases in magnitude eventually canceling the scalar en¬ 
ergy precipitating a bounce. After the bounce as the 
universe starts to inhate the curvature term vanishes, see 
Fig© 

B. Perturbations in a closed universe 

If the pre-bounce dehationary phase is allowed to con¬ 
tinue all the way to past inhnity, then all the perturbative 



0 1 2 3 4 5 


FIG. 6: Background Evolution in Bounce-inflation model. 
The black, red and blue curves represent a(t), pc/>(t) and the 
absolute value of the curvature density, A/ 4 /a 2 (t), respec¬ 
tively. 


modes would essentially start out superHubble. At some 
point during contraction they all become subHubble ex¬ 
cept the really long wavelength modes which remain 
longer than the bounce length scale even at the bounce. 
After the bounce, exactly the reverse happens, the modes 
which became subHubble starts to exit outside the Hub¬ 
ble radius as the universe inhates. In such a scenario it 
is not clear what initial conditions one should choose for 
the various fluctuations as they start in the superHub¬ 
ble phase! For instance, the all-important inhationary 
prediction of a near scale-invariant spectrum of pertur¬ 
bations that relies on being able to impose the Bunch 
Davis (BD) vacuum initial conditions in the subHubble 
phase, would no longer be applicable. Accordingly, phe¬ 
nomenologically a more interesting scenario would be to 
have a non-dehationary (say, matter or radiation domi¬ 
nated) contraction phase prior to dehation. In this case 
all modes would start out subHubble where the BD ini¬ 
tial conditions can be imposed. Some of them will exit 
the Hubble radius during this phase. Once the deha¬ 
tionary phase begins though, all the modes that haven’t 
exited will continue their journey as subHubble modes 
all the way to the expansion phase. It is possible there¬ 
fore for these subHubble modes to exit the Hubble radius 
only once during the familiar inhationary phase, thereby 
reproducing the near-scale-invariance of the ensuing pri¬ 
mordial spectrum. Can this process work? 

Among the modes that do become superHubble dur¬ 
ing contraction, some of them re-enter the Hubble hori¬ 
zon during the dehationary phase, while others remain 
superHubble forever. What kind of transformations do 
these perturbations undergo during the bounce? Clearly 
these are key questions for various model building exer¬ 
cises. Our aim in this paper is not to try and construct 
models, but rather to shed light on how perturbations 
(sub and super) evolve during the bounce which will aid 
for future model building. As we will see, our study will 
also provide us ways to go beyond the formalism and 
results we have discussed in the earlier sections. 
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In dealing with fluctuations from scalar fields, often 
it is convenient to work with the canonical variable \ 
defined via (IIII.2dl) where 


n 


6d> 

$ + 77^2 


z 

K c 


and ei 


—ay/ 2ei[l + K c e\/(k 2 
M 4 



2 H 2 M% 2 H 2 M 2 ' 


4 K c )]~ 4 


(VI.67) 


Above, 5(j) refers to the fluctuation of the scalar field. 

In [693 the evolution equation of \ was derived for a 
closed universe: 


Xk + 


k 2 - 3 K r - 2 


K c z' 

Hz 


z 

z 


Xk = o (VI.68) 


We observe is that at the bounce point when 71 = 0, 
the fourth term in the evolution equation blows up. This 
is sometimes regarded as evidence of the fact that per¬ 
turbations grow in an unbounded manner at the bounce 
leading one to argue that bouncing cosmologies are phys¬ 
ically untenable or problematic. On the contrary, what 
this really indicates is simply that x is n °t a good variable 
to use at the bounce point! Instead if one uses, <h, 8<j>, to 
track the fluctuations one finds completely smooth tran¬ 
sitions around the bounce, just as we saw around the 
turnaround cosmology in section IIII Bl 

To demonstrate this, let us look at the evolution equa¬ 
tions for <f>, a = M p 8(j)/(j)' in a closed universe [69j |: 



^' k + n^ k 


a' k + Ha k 


4 > 


2 Ml 


1 + 


2M 2 (—k 2 + 3 K c 


(VI.69) 
<!>,, (VITO) 


FIG. 7: Plots of $ (top), a (middle) and \ (bottom) for 
different comoving wavenumbers, k = 0.01 (red), 1 (green), 
10 (blue) and 100 (purple) for the case M = M p = 1. Our 
initial conditions are a(0) = 0.1, 0(0) = 1 and (p'(0) = 0. 


By inspection it is clear that the above equations remain 
regular during the bounce and can indeed be solved to 
yield perfectly regular solutions for <h, a. Fig. [HI top, 
middle and bottom represent the evolution of <f>, a and x 
respectively. We notice that all the variables, including 
X are finite and well defined. However x is discontinuous 
at the bounce, and this is the reason why it cannot be 
tracked directly from its evolution equation (IVI.68I) . It’s 
also obvious now that this does not represent any phys¬ 
ical instability and one just has to use the appropriate 
perturbative variable to track the fluctuations across the 
bounce. In our opinion, this is quite an important result 
and along with our similar findings for the turnarounds 
should go a long way to put to rest some of the concerns 
regarding the consistency of bouncing cosmologies. Let 
us now try to glean out some general patterns of evolu¬ 
tion of both the sub and super Hubble fluctuations. 


C. Evolution of Fluctuations 

1. SubHubble Evolution 

In Fig. [8] we have plotted <h, a and x as functions of 
time for different comoving fluctuations. As is evident 
from inspection, when the wavelength becomes smaller 
than the bounce scale (i.e. subHubble) the amplitude 
of oscillations for the x variable freeze. This is exactly 
what was found for the turnaround phase as well, see 
section UlI Bl Indeed such a behavior is well known from 
analysis of subHubble perturbations in inflationary cos¬ 
mology, but it is rather pleasing to see that the trend 
continues in more complex cosmologies lending credence 
to the fact that this may be a general feature of subHub¬ 
ble fluctuations. 
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What can one expect if and when one includes quan¬ 
tum gravitational effects that should be relevant for small 
wave-length fluctuations? SubHubble perturbations in 
FLRW backgrounds are actually very similar to fluctua¬ 
tions in Minkowski space-time as the wavelengths are too 
small to see the cosmological expansion. The only effect 
of the expansion is to stretch the wavelengths. This corre¬ 
spondence becomes particularly visible in conformal time 
when one works with the “canonical variable”, X- Barring 
near the bounce, it simply undergoes sinusoidal oscilla¬ 
tions. In ii a very general class of higher derivative 
theories of gravity, which may be able to incorporate 
several QG effects, have been studied. In particular it 
was shown that the perturbations break up into different 
modes, each mode effectively evolving according to its 
own Klein-Gordon type equation with a particular mass. 
This strongly suggests that in the subHubble regime we 
may continue to see simple sinusoidal oscillatory behavior 
in the appropriate canonical variables. In other words, 
we expect to only see an adiabatic stretching of the wave¬ 
lengths of the subHubble fluctuations with the evolution 
of the universe, while its amplitude remains a constant. 

This would, for instance, strongly support the idea that 
even if the standard inflationary phase is preceded by 
a bounce, the l/Vk dependence of the “initial” ampli¬ 
tude of subHubble fluctuations, which comes from im¬ 
posing the Bunch-Davis initial conditions in the contract¬ 
ing phase, will be preserved in the course of the bounce, 
thereby enabling the inflationary mechanism to provide 
the observed near scale-invariant seed fluctuations for 
CMB, as usual. 


2. SuperHubble Evolution 

What happens to the superHubble fluctuations? The 
above model actually represents a case where neither of 
our superHubble formulas, (III. 231) and (III. 301) . are appli¬ 
cable. (III.231) is not applicable because we have a closed 
universe where the scale-factor symmetry is no longer 
valid, while (III. 301) works for a closed universe but only 
for ideal fluids; we have a scalar field. In fact, we espe¬ 
cially chose this model to, on one hand, show the limi¬ 
tations of our superHubble formalism, but on the other, 
suggest how progress can still be made in spite of it. For 
such cases, one has to know the perturbative evolutions 
equations which can be solved numerically. What is use¬ 
ful though is that our formalism can still help in obtaining 
the transfer matrix that connects the modes in contrac¬ 
tion phase with modes during expansion. This is because 
we know that even in this case, once the curvature term 
becomes small, our superHubble formalism becomes ap¬ 
plicable and therefore the perturbations should still be¬ 
have according to (IIV.30D at late times. We thus expect 
that even in these models where our approach cannot 
directly be used to track the superHubble fluctuations 
during the entire cosmological evolution, with the help 
of numerics, we can still get a handle on the evolution 


of these fluctuations. Let us see this explicitly in our 
specific model now. 

We know that at large times in the past and future, 
the spatial curvature is negligible and therefore the per¬ 
turbation equation for $ must reduce to the standard 
form I7l|: 

$fc+2(? 7 -e)'H$ , fc + fc 2 $ fc + 2(77-2e)'H 2 $ fc = 0, (VI.71) 


where the scalar field is assumed to be “slowly rolling”, 
and 


e = 


1 = 



(VI.72) 
(VI.73) 


are the slow-roll parameters. For our simple potential 
(IVL641) . as long as 


M, 


Vl 


(VI.74) 


the slow-roll approximation should be valid with 


e ~ and ?? = 0 


2u 2 


(VI.75) 


For an exponential contraction or expansion the scale 
factor is given by 


a(r) = =F 


1 


A(r - T±) 


and EL = -f - 


1 


T — T± 


(VI.76) 


where r —> r+ from below for expansion and r — > t_ from 
above for contraction. In either case, the perturbation 
equation can be solved in the superHubble regime to yield 


$ = C ± (M p t)‘* +£>±(M p r)« 


where, 


91 

92 


4e 

1 - 2e 


0 


1 - 2e + 


4e 

1 - 2e 


1 


(VI. 77) 


(VI. 78) 
(VI. 79) 


In other words, approximately one can rewrite the super¬ 
Hubble modes as 

$ = c ± + , (VI.80) 

a(r) 

precisely the form one obtains from our analytical for¬ 
mula (IIV.30I) . 

We can check this explicitly by solving the exact dif¬ 
ferential equations (IVI.70I) and then matching it with 
(IVI.80I) for large times. To see this asymptotic match¬ 
ing, it is more convenient to work with proper time as 
opposed to conformal time (the universe ends and begins 
abruptly in conformal time). In Fig. [5J we have plotted 



























18 



FIG. 8: The blue curve represents 4>(a) that we obtain nu¬ 
merically. The yellow and green curves represent the analytic 
approximations in the far past and distant future respectively. 
Our initial conditions are a(0) = 0.1, </>(0) = 1 and = 0. 


the evolution of $ as the universe undergoes a bounce. 
As one can see, at large times, both in the past and fu¬ 
ture, it is fitted extremely well with (1VI.80I) , represented 
by green and yellow curves respectively. 

While fitting the asymptotic forms we could easily ob¬ 
tain the coefficients C ±, D ±, and we noticed a somewhat 
surprising result. It seems that the constant mode ampli¬ 
tude actually changes during the course of the bounce as 
evident in Fig. [8] as well. Also, the change in amplitude 
seems to vary with k\ This means that the spectrum of 
the constant mode could be getting altered by contribu¬ 
tions from the growing mode of the contraction phase. 
This is very different from the flat universe case where 
no such alteration was possible. This may be good news 
for models based on ekpyrotic cosmology which natu¬ 
rally tend to produce a scale-invariant spectrum in the 
growing mode, and warrants a more comprehensive in¬ 
vestigation in the future. Our analysis can explains is 
why at times different results were obtained in different 
models when addressing the question of mode mixing: A 
first analysis of the evolution of fluctuations through a 
specific nonsingular bounce obtained by using a higher 
derivative gravity action was made in [471 ] . showing that 
initial scale-invariant fluctuations do not pass through 
the bounce, thus confirming the results of [zlHzli]- On 
the other hand, bounces obtained by adding spatial cur¬ 
vature and matter with wrong-sign kinetic terms, but 
with the standard gravitational action, have been studied 
in ErElzl, with differing results. It would also be most 
interesting to see how our analysis and results gets ex¬ 
tended/changed when one ventures to more complicated 
cosmology involving many different degrees of freedom! 


VII. CONCLUSIONS 

Let us conclude by highlighting the main results of our 
paper and suggesting a few promising future research di¬ 


rections. We have established a comprehensive connec¬ 
tion between the evolution of the superHubble fluctua¬ 
tions and background cosmology. This is manifested 

• in our ability to obtain differential equations for 
the evolution of superHubble fluctuations from the 
evolution equations of background quantities like 
the scale-factor and various energy densities, as il¬ 
lustrated in the appendix 

• in our ability to directly obtain the superHubble 
modes from the knowledge of all possible back¬ 
ground solutions, and 

• in our being able to derive explicit analytical for¬ 
mulas for two of the usual superHubble modes as¬ 
sociated with the scale-factor and time-translation 
symmetry. 

We were able to verify our results by explicitly solving 
perturbation equations analytically/numerically in some 
special toy models in GR where the perturbation equa¬ 
tions were known. We not only looked at monotonically 
expanding and contracting universes but also more com¬ 
plex cosmologies involving bounces and turnarounds me¬ 
diating nonsingular transitions between the contracting 
and expanding phases. This is particularly relevant for 
bouncing and cyclic cosmologies which were the main fo¬ 
cus of the paper. 

In the process, 

• we verified that perturbations can indeed remain 
small during the course of bounces and turnarounds 
as long as one is careful to use appropriate vari¬ 
ables, divergences in ill-defined variables can indeed 
be misleading. 

• We were able to prove that for a spatially flat uni¬ 
verse, the constant mode in the expanding branch 
exclusively couples to the constant mode in the con¬ 
tracting branch confirming some of the earlier re¬ 
sults. However in a closed universe the transfer 
matrix may be more interesting. This latter result 
needs to be investigated more comprehensively in 
the future due to it’s obvious phenomenological im¬ 
plications. 

• Our numerical forays allowed us to corrobo¬ 
rate the intuitive argument j53j, !6Cj that al¬ 
though technically all modes are subHubble at the 
bounce/turnaround points, they behave as super¬ 
Hubble modes as long as their wavelengths are 
larger than the appropriate cosmological time scale 
controlled essentially by the energy density scale of 
the bounces and turnarounds in question. 

• During our numerical analysis we found that the 
amplitude of the sinusoidal oscillations in the 
canonical perturbation variable, y, remains a con¬ 
stant for the subHubble modes throughout the cos¬ 
mological evolution. The fact that they should re¬ 
main a constant in monotonically expanding and 
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contracting phases was well known, but that it 
does so around the bounce and the turnaround is 
a pleasant surprise, especially considering that x is 
ill-defined precisely at the bounce and turnaround. 
We provided further arguments to suggest why this 
may be a very generic feature even when quantum 
gravitational effects are introduced. 

• After testing our formalism in simple toy models, 
we applied our results to the emergent cyclic infla¬ 
tion scenario. In brief, our results indicated that 
superHubble fluctuations grow in emergent cyclic 
phase, so that one can only envision homogenous 
patches of Planck length. This homogeneity is how¬ 
ever preserved in the cyclic inflation phase conform¬ 
ing the viability of this model. 

• Next we applied our formalism to a bounce inflation 
model. We first checked that perturbations can re¬ 
main finite and small during the bounce. Moreover, 
we found that the subHubble fluctuations can pre¬ 
serve its initial Bunch-Davis spectrum so that the 
generation of the near-scale-invariant spectrum of 
fluctuations in inflationary paradigm can operate 
in the usual way. 

Throughout the course of this paper we have seen how 
our work and results can be generalized in different di¬ 
rections. Perhaps the most urgent and interesting ven¬ 
ture would be to study a multi-component cosmological 
system and in particular study the possible generation 
of near scale-invariant perturbations in cosmologies such 
as new ekpyrotic, matter-bounce and lragedorn-bounce 
models. Another important direction would be to try 
to come up with a general framework to understand the 
subHubble fluctuations which, together with our results 
on superHubble fluctuations, would then provide a very 
concrete basis to construct and validate/reject bounc¬ 
ing/cyclic cosmological paradigms. 

VIII. APPENDIX 

In this appendix we illustrate how one may be able to 
derive the evolution equations that superHubble modes 
must satisfy from the knowledge of the background evo¬ 
lution equations by looking at the simple example of a 
single ideal fluid in GR. So, our starting point are the 
Eqns. (IIII. 171) and (III.271) describing the evolution of 
the scale-factor and energy density in a homogeneous 
isotropic cosmological background. Now, small tempo¬ 
ral fluctuations of the scale factor and the energy density 
is captured by 

U = U + ip' (VIII.81) 

and p = p(l + 5) . (VIII.82) 

Then, perturbing the Hubble equation one obtains 

2 USU = K[2ap(5a) + a 2 {5p)} 


=> 2ip' = U{2 i/j + 5) . (VIII.83) 

Similarly, perturbing the continuity equation one obtains 
(Sp)' + 3(1 + co)[p(SU) + U(Sp)} 


— S\p' + 3(1 + uj)Up\ + p\S' + 3(1 + w)(<7H)] — 0 

where we have used the background continuity equation 
(III.271) . Finally we have 

S' + 3ip'( 1 + w) = 0 (VIII.84) 

(IVITT.S3D and (IVTTT.S4I) can now be combined to obtain 
a single differential equation for if>: 

i/j" + U(1 + 3u)ip' = 0 (VIII.85) 

We can now substitute '0 in terms of and £. For 
instance, we find 

ip' = - 2- £H' 

ip" = - 2&U - -mu' - iU" (VIII.86) 

Then using (III. 91) . and (IVIII.86I) We find 

$"+3(l+uj)H<S> , +2[2H'+(l+3uj)H 2 }$+Z['H"+nH'(l+3Lo)\ = 0 

which reduces to the super-Hubble limit (IIII.lHl) using 
the identity 

2 U' + (1 + 3 u>)H 2 = 0 (VIII.87) 

valid in GR. 

This demonstrates how starting from the homoge¬ 
neous isotropic equations one can deduce the perturba¬ 
tion equation obeyed by $ in the super-Hubble phase. 

A similar algorithm can be followed to find superHub¬ 
ble fluctuations from the background equations for any 
covariant theory of gravity. Note, that this procedure 
does not rely on symmetries such as the time translation 
or scale-factor symmetry and is therefore more general 
than the approach we have mostly followed in our paper. 

It should also be possible to generalize the algorithm to 
include scalar fields but we leave this as a future exercise. 

Acknowledgments: TB would like to thank Abhay 
Ashtekar for initial collaboration and several helpful dis¬ 
cussions especially on the subject of Loop Quantum cor¬ 
rections to cosmological equations. TB would also like 
to thank Anupam Mazumdar and Robert Brandenberger 
for several insightful discussions on cosmological pertur¬ 
bation theory. 


















20 


[1] S. Hawking and R. Penrose, “The Singularities of grav¬ 
itational collapse and cosmology,” Proc.Roy.Soc.Lond., 
vol. A314, pp. 529-548, 1970. 

[2] A. Borde and A. Vilenkin, “Eternal inflation and the ini¬ 
tial singularity,” Phys.Rev.Lett., vol. 72, pp. 3305-3309, 
1994. 

[3] A. Borde, A. H. Guth, and A. Vilenkin, “Infla¬ 
tionary space-times are incompletein past directions,” 
Phys.Rev.Lett., vol. 90, p. 151301, 2003. 

[4] G. F. Ellis and R. Maartens, “The emergent uni¬ 
verse: Inflationary cosmology with no singularity,” 

Class. Quant. Grav., vol. 21, pp. 223-232, 2004. 

[5] G. F. Ellis, J. Murugan, and C. G. Tsagas, 
“The Emergent universe: An Explicit construction,” 
Class.Quant.Grav., vol. 21, pp. 233-250, 2004. 

[61 T. Biswas, “The Hagedorn Soup and an Emergent Cyclic 
Universe,” 2008. 

[7] A. Conroy, A. S. Koshelev, and A. Mazumdar, “Geodesic 
completeness and homogeneity condition for cosmic infla¬ 
tion,” Phys.Rev., vol. D90, no. 12, p. 123525, 2014. 

[8] M. Gasperini and G. Veneziano, “The Pre - big bang sce¬ 
nario in string cosmology,” Phys.Rept., vol. 373, pp. 1- 
212, 2003. 

[9] J. Khoury, B. A. Ovrut, P. J. Steinhardt, and N. Turok, 
“The Ekpyrotic universe: Colliding branes and the origin 
of the hot big bang,” Phys.Rev., vol. D64, p. 123522, 
2001 . 

[10] F. Finelli and R. Brandenberger, “On the generation of a 
scale invariant spectrum of adiabatic fluctuations in cos¬ 
mological models with a contracting phase,” Phys.Rev., 
vol. D65, p. 103522, 2002. 

[11] R. H. Brandenberger, “The Matter Bounce Alternative 
to Inflationary Cosmology,” 2012. 

[12] A. Nayeri, R. H. Brandenberger, and C. Vafa, “Produc¬ 
ing a scale-invariant spectrum of perturbations in a Hage¬ 
dorn phase of string cosmology,” Phys.Rev.Lett., vol. 97, 
p. 021302, 2006. 

[13] T. Biswas, R. Brandenberger, A. Mazumdar, and 
W. Siegel, “Non-perturbative Gravity, Hagedorn Bounce 
& CMB,” JCAP, vol. 0712, p. Oil, 2007. 

[14] T. Biswas, T. Koivisto, and A. Mazumdar, “Atick- 
Witten Hagedorn Conjecture, near scale-invariant matter 
and blue-tilted gravity power spectrum,” 2014. 

[15] T. Biswas and A. Mazumdar, “Super-Inflation, 
Non-Singular Bounce, and Low Multipoles,” 
Class.Quant.Grav., vol. 31, p. 025019, 2014. 

[16] Z.-G. Liu, Z.-K. Guo, and Y.-S. Piao, “Obtaining the 
CMB anomalies with a bounce from the contracting 
phase to inflation,” Phys.Rev., vol. D88, p. 063539, 2013. 

[17] G. Calcagni, L. Modesto, and P. Nicolini, “Super- 
accelerating bouncing cosmology in asymptotically-free 
non-local gravity,” Eur.Phys.J., vol. C74, no. 8, p. 2999, 

2014. 

[18] Y.-T. Wang and Y.-S. Piao, “Parity violation in pre¬ 
inflationary bounce,” Phys.Lett., vol. B741, pp. 55-60, 

2015. 

[19] Y. Cai, Y.-T. Wang, and Y.-S. Piao, “Pre-inflationary 
primordial perturbations,” 2015. 

[20] C. Bennett et al, “Nine-Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Final Maps 
and Results,” Astrophys.J.Suppl., vol. 208, p. 20, 2013. 


[21] P. Ade et al., “Planck 2013 results. XVI. Cosmological 
parameters,” 2013. 

[22] R. Brandenberger, T. Biswas, T. Koivisto, and 
A. Mazumdar, “Cosmological perturbations from Statis¬ 
tical Fluctuations,” 2013. 

[23] Y.-F. Cai, C. Gao, and E. N. Saridakis, “Bounce and 
cyclic cosmology in extended nonlinear massive gravity,” 
JCAP, vol. 1210, p. 048, 2012. 

[24] T. Prokopec, “A Solution to the cosmological constant 
problem,” 2006. 

[25] Y.-S. Piao, “Design of a Cyclic Multiverse,” Phys.Lett., 
vol. B691, pp. 225-229, 2010. 

[26] Y.-S. Piao, “Proliferation in Cycle,” Phys.Lett., 
vol. B677, pp. 1-5, 2009. 

[27] Y.-S. Piao, “Can the universe experience many cycles 
with different vacua?,” Phys.Rev., vol. D70, p. 101302, 
2004. 

[28] D. J. Mulryne, N. Nunes, R. Tavakol, and J. E. Lid- 
sey, “Inflationary cosmology and oscillating universes in 
loop quantum cosmology,” Int.J.Mod.Phys., vol. A20, 
pp. 2347-2357, 2005. 

[29] A. Ashtekar, T. Pawlowski, and P. Singh, “Quantum Na¬ 
ture of the Big Bang: Improved dynamics,” Phys.Rev., 
vol. D74, p. 084003, 2006. 

[301 M. Boiowald, “Loop quantum cosmology,” Livinq 
Rev.Rel, vol. 8, p. 11, 2005. 

[31] Y. Shtanov and V. Sahni, “Bouncing brane worlds,” 
Phys.Lett., vol. B557, pp. 1-6, 2003. 

[32] K. Freese, M. G. Brown, and W. H. Kinney, “The Phan¬ 
tom Bounce: A New Proposal for an Oscillating Cosmol¬ 
ogy,” 2008. 

[33] L. Baum and P. H. Frampton, “Turnaround in cyclic cos¬ 
mology,” Phys.Rev.Lett., vol. 98, p. 071301, 2007. Altered 
and improved model astro-ph/0608138 with new title. 

[34] T. Biswas, A. Mazumdar, and W. Siegel, “Bouncing 
universes in string-inspired gravity,” JCAP, vol. 0603, 
p. 009, 2006. 

[35] T. Biswas, E. Gerwick, T. Koivisto, and A. Mazumdar, 
“Towards singularity and ghost free theories of gravity,” 
Phys.Rev.Lett., vol. 108, p. 031101, 2012. 

[36] T. S. Koivisto, “Bouncing Palatini cosmologies and their 
perturbations,” Phys.Rev., vol. D82, p. 044022, 2010. 

[37] Y.-F. Cai and E. N. Saridakis, “Non-singular Cyclic Cos¬ 
mology without Phantom Menace,” J.Cosmol., vol. 17, 
pp. 7238-7254, 2011. 

[38] W. Israel, “Singular hypersurfaces and thin shells in gen¬ 
eral relativity,” Nuovo Cim., vol. B44S10, p. 1, 1966. 

[39] J.-c. Hwang and E. T. Vishniac, “Gauge-invariant joining 
conditions for cosmological perturbations,” Astrophys.J., 
vol. 382, pp. 363-368, 1991. 

[40] N. Dcruelle and V. F. Mukhanov, “On matching 
conditions for cosmological perturbations,” Phys.Rev., 
vol. D52, pp. 5549-5555, 1995. 

[41] J. Martin, P. Peter, N. Pinto Neto, and D. J. Schwarz, 
“Passing through the bounce in the ekpyrotic models,” 
Phys.Rev., vol. D65, p. 123513, 2002. 

[42] P. Peter and N. Pinto-Neto, “Primordial perturbations 
in a non singular bouncing universe model,” Phys.Rev., 
vol. D66, p. 063509, 2002. 

[43] P. Peter, N. Pinto-Neto, and D. A. Gonzalez, “Adia¬ 
batic and entropy perturbations propagation in a bounc- 



21 


ing universe,” JCAP, vol. 0312, p. 003, 2003. 

[44] J. Martin and P. Peter, “On the causality argu¬ 
ment in bouncing cosmologies,” Phys.Rev.Lett., vol. 92, 
p. 061301, 2004. 

[45] J. Martin and P. Peter, “On the properties of the 
transition matrix in bouncing cosmologies,” Phys.Rev., 
vol. D69, p. 107301, 2004. 

[46] R. Durrer and F. Vernizzi, “Adiabatic perturbations in 
pre - big bang models: Matching conditions and scale 
invariance,” Phys.Rev., vol. D66, p. 083503, 2002. 

[47] S. Tsujikawa, R. Brandenberger, and F. Finelli, “On the 
construction of nonsingular pre - big bang and ekpyrotic 
cosmologies and the resulting density perturbations,” 
Phys.Rev., vol. D66, p. 083513, 2002. 

[48] D. Salopek and J. Bond, “Nonlinear evolution of long 
wavelength metric fluctuations in inflationary models,” 
Phys.Rev., vol. D42, pp. 3936-3962, 1990. 

[49] M. Sasaki and T. Tanaka, “Superhorizon scale dynam¬ 
ics of multiscalar inflation,” Prog. Theor.Phys., vol. 99, 
pp. 763-782, 1998. 

[50] H. Kodama and T. Hamazaki, “Evolution of cosmological 
perturbations in the long wavelength limit,” Phys.Rev., 
vol. D57, pp. 7177-7185, 1998. 

[51] T. Biswas and S. Alexander, “Cyclic Inflation,” 
Phys.Rev., vol. D80, p. 043511, 2009. 

[52] T. Biswas and A. Mazumdar, “Inflation with a negative 
cosmological constant,” Phys.Rev., vol. D80, p. 023519, 
2009. 

[53] T. Biswas, A. Mazumdar, and A. Shafieloo, “Wiggles 
in the cosmic microwave background radiation: echoes 
from non-singular cyclic-inflation,” Phys.Rev., vol. D82, 
p. 123517, 2010. 

[54] T. Biswas, T. Koivisto, and A. Mazumdar, “Could our 
Universe have begun with Negative Lambda?,” 2011. 

[55] W. Duhe and T. Biswas, “Numerically Investigating the 
Emergent Cyclic Inflation Scenario,” Class. Quant. Grav., 
vol. 31, p. 155010, 2014. 

[56] T. Biswas, “Before the Bang,” arXiv[gr-qc] 1308.0024, 
2013. 

[57] A. Cardoso and D. Wands, “Generalised perturbation 
equations in bouncing cosmologies,” Phys.Rev., vol. D77, 
p. 123538, 2008. 

[58] V. F. Mukhanov, H. Feldman, and R. H. Brandenberger, 
“Theory of cosmological perturbations. Part 1. Classical 
perturbations. Part 2. Quantum theory of perturbations. 
Part 3. Extensions,” Phys.Rept., vol. 215, pp. 203-333, 
1992. 

[59] T. Biswas, T. Koivisto, and A. Mazumdar, “Phase transi¬ 
tions during Cyclic-Inflation and Non-gaussianity,” 2013. 

[60] G. Geshnizjani, W. H. Kinney, and A. M. Dizgah, “Gen¬ 
eral conditions for scale-invariant perturbations in an ex¬ 
panding universe,” JCAP, vol. 1111, p. 049, 2011. 


[61] T. Biswas, R. Brandenberger, D. A. Easson, and 
A. Mazumdar, “Coupled inflation and brane gases,” 
Phys.Rev., vol. D71, p. 083514, 2005. 

[62] L. Hui and W. H. Kinney, “Short distance physics and 
the consistency relation for scalar and tensor fluctua¬ 
tions in the inflationary universe,” Phys.Rev., vol. D65, 
p. 103507, 2002. 

[63] J. Magueijo and P. Singh, “Thermal fluctuations in loop 
cosmology,” Phys.Rev., vol. D76, p. 023510, 2007. 

[64] Y.-F. Cai, W. Xue, R. Brandenberger, and X.-m. 
Zhang, “Thermal Fluctuations and Bouncing Cosmolo¬ 
gies,” JCAP, vol. 0906, p. 037, 2009. 

[65] C. Hull, “De Sitter space in supergravity and M theory,” 
JHEP , vol. 0111, p. 012, 2001. 

[66] S. Kachru et al ., “Towards inflation in string theory,” 
JCAP, vol. 0310, p. 013, 2003. 

[67] M. R. Douglas and S. Kachru, “Flux compactification,” 
Rev. Mod. Phys., vol. 79, pp. 733-796, 2007. 

[68] A. Mazumdar and J. Rocher, “Particle physics models of 
inflation and curvaton scenarios,” Phys.Rept., vol. 497, 
pp. 85-215, 2011. 

[69] D.-H. Zhang and C.-Y. Sun, “The Exact evolution equa¬ 
tion of the curvature perturbation for closed universe,” 
Chin.Phys.Lett., vol. 21, p. 1865, 2004. 

[70] T. Biswas, A. S. Koshelev, A. Mazumdar, and S. Y. 
Vernov, “Stable bounce and inflation in non-local higher 
derivative cosmology,” JCAP, vol. 1208, p. 024, 2012. 

[71] A. Riotto, “Inflation and the theory of cosmological per¬ 
turbations,” pp. 317-413, 2002. 

[72] D. H. Lyth, “The Primordial curvature perturbation in 
the ekpyrotic universe,” Phys.Lett., vol. B524, pp. 1-4, 
2002 . 

[73] R. Brandenberger and F. Finelli, “On the spectrum of 
fluctuations in an effective field theory of the Ekpyrotic 
universe,” JHEP, vol. 0111, p. 056, 2001. 

[74] J.-c. Hwang, “Cosmological structure problem in the 
ekpyrotic scenario,” Phys.Rev., vol. D65, p. 063514, 2002. 

[75] J. Khoury, B. A. Ovrut, N. Seiberg, P. J. Steinhardt, 
and N. Turok, “From big crunch to big bang,” Phys.Rev., 
vol. D65, p. 086007, 2002. 

[76] P. Creminelli, A. Nicolis, and M. Zaldarriaga, “Pertur¬ 
bations in bouncing cosmologies: Dynamical attractor 
versus scale invariance,” Phys.Rev., vol. D71, p. 063505, 
2005. 

[77] F. Finelli, “Study of a class of four-dimensional nonsin¬ 
gular cosmological bounces,” JCAP, vol. 0310, p. Oil, 
2003. 

[78] L. E. Allen and D. Wands, “Cosmological perturba¬ 
tions through a simple bounce,” Phys.Rev., vol. D70, 
p. 063515, 2004. 



